Tunable Excitons in Biased Bilayer Graphene 
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Recent measurements have shown that a continuously tunable bandgap of up to 250 meV can 
be generated in biased bilayer graphene [Y. Zhang et at., Nature 459, 820 (2009)], opening up 
pathway for possible graphene-based nanoelectronic and nanophotonic devices operating at room 
temperature. Here, we show that the optical response of this system is dominated by bound excitons. 
The main feature of the optical absorbance spectrum is determined by a single symmetric peak 
arising from excitons, a profile that is markedly different from that of an interband transition picture. 
Under laboratory conditions, the binding energy of the excitons may be tuned with the external 
bias going from zero to several tens of meV's. These novel strong excitonic behaviors result from a 
peculiar, effective "one-dimensional" joint density of states and a continuously-tunable bandgap in 
biased bilayer graphene. Moreover, we show that the electronic structure (level degeneracy, optical 
selection rules, etc.) of the bound excitons in a biased bilayer graphene is markedly different from 
that of a two-dimensional hydrogen atom because of the pseudospin physics. 
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The low-energy electronic states of graphene are de- 
scribed by a massless Dirac equation [lH3|. If an ex- 
tra layer is added [Fig. [B^a)], the electronic properties 
change drastically and the charge carriers become mas- 
sive [Fig.[ljb)] [45. There have been a number of theoret- 
ical studies on the possibility of opening up a bandgap 
in the gapless bilayer graphene if an electric field is ap- 
plied perpendicularly [Figs.[Uc) and 1(d)] [B-15|. Indeed, 
a bandgap has been observed in the case of an internal 
perpendicular electric field generated by an imbalance of 
doped charge between the two graphene layers jl6|. Also, 
bandgap opening in bilayer graphene under an electric 
field from a sin g l e g ate has been observed by infra-red 



spectroscopy 17h21|. quantum Hall measurement 



and scanning tunneling spectroscopy [23| . 

A bandgap opening up in bilayer graphene under an 
electric field from a double-gate configuration has fur- 
ther been observed in transport experiments [l^, [IH] and 
quantum Hall measurements [26| . Very recently, infra- 
red measurements showed that the bandgap of bilayer 
graphene in a double-gate geometry is continuously tun- 
able up to 250 meV, an order of magnitude higher than 
the thermal energy at room temperature (27| . This dis- 
covery provides exciting new possibilities for the nano- 
electronic and nanophotonic device applications of bi- 
layer graphene at room temperature. 

Theoretical studies on the optical response of intrinsic 



bilayer graphene within a single-particle picture [28|, 29 1 
as well as including electron- hole {e-h) interactions 30, 
3l| have been performed. It is found that there are neg- 
ligible many-electron effects on the low-energy (< 1 eV) 
optical response of graphene and bilayer graphene [30| . 
There have also been theoretical studies within a single- 
particle picture on the electronic and optical properties 
of biased bilayer graphene (BBG) ,32-35]. However, the- 
oretical investigation of excitonic effects on the optical 
response of this novel tunable bandgap system has yet to 



be performed up to now. It is known that e-h interactions 
play a crucial role in the optical response of semiconduc- 
tors [IB], especially, semiconducting nanostructures [3?! - 
|40| . Excitonic effects in BBG with a finite bandgap are 
expected to be important, considering that the lowest- 
energy van Hove singularity in its joint electronic density 
of states exhibits a one-dimensional (ID), and not a two- 
dimensional (2D), behavior (i. e. , it diverges as inverse of 
the energy difference from the bandgap) .32] . 

Here, we obtain the optical response of a BBG includ- 
ing e-h interactions by solving the Bethe-Salpeter equa- 
tion (BSE): 
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where A^^^^ is the amplitude of a free e-h pair configu- 
ration composed of the electron state ]ck) and the hole 
state Iwk), fl^ is the exciton excitation energy, Eck and 
Ey]i are quasiparticle energies, and K'^^ is the e-h inter- 
action kernel [36| . The absorption spectrum is calculated 
by evaluating the optical matrix elements [sil ] using the 
eigenstates and eigenvalues of the BSE. 

As in recent experiments 27|, |4l|, we focus here on 
the case in which the net charge on the BBG is zero, or, 
the displacement fields D above and below the bilayer 
graphene are the same [Fig. [He)] . We find that the opti- 
cal response of BBG is dominated by low-energy bound 
excitons with huge oscillator strength due to the ID na- 
ture in the joint density of states. As a consequence, the 
main peak of the absorbance profile becomes highly sym- 
metric. The binding energy and oscillator strength of the 
excitons increase with the bandgap. We find a very rich 
electronic structure for the excitons in a BBG. Especially, 
we discover a symmetry breaking of excitons having an- 
gular momenta of equal magnitude but of opposite sign 
which leads to an unusual selection rule in the optical 
absorption. This phenomenon is explained in terms of 
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the pseudospin, a degree of freedom describing the bond- 
ing character between neighboring carbon atoms 3] , in a 
BBG. 

In this study, we make use of the k -p based method de- 
veloped by Ando and coworkers for the excitonic spectra 
of graphene and carbon nanotubes j42l446l | . Although, 
unlike the first-principles GM^-BSE approach [s^ that is 
parameter free, the current method is based on a tight- 
binding formalism and treats electron-electron interac- 
tions within the screened Hartree-Fock approximation, 
it does provide excitonic features of the absorption pro- 
file that may be compared with experiments for complex 
structures and applied fields 4J, |47[. For the e-h kernel 
K'^^, we consider only the attractive direct term, which 
is by far dominant and describes the screened interaction 
between electrons and holes, and neglect the repulsive 
exchange term. The exchange kernel is responsible for 
singlet-triplet splitting and the splitting among states 
within individual singlet and triplet complexes, but is 
usually only a few percent in magnitude of the direct 
term 

The quasiparticle energies i?ck = £ck + Sck and i?„k — 
e„k + S„k are obtained by first calculating the bare en- 
ergy £k within the k ■ p formalism ji^ l using a tight- 
binding Hamiltonian where we set the intralayer hop- 
ping parameter between the nearest-neighboring atoms 

70 = 2.6 eV and the interlayer hopping parameter 

71 — 0.37 eV. These parameters reproduce well the 
bandstructure of pristine bilayer graphene obtained from 
density-functional calculations within the local density 
approximation (LDA) [s^]- The self energy Ek is calcu- 
lated within the screened Hartree-Fock approximation, 
using the static random-phase dielectric function [42- 
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47[. We calculate the static polarizability within 
the random-phase approximation 42 - 3 by includ- 



ing the four electronic bands closest to the bandgap aris- 
ing from the tt states with an energy cutoff of 5 eV 
(we have checked that the resulting quasiparticle en- 
ergies are insensitive to this cutoff), and incorporate 
the effects of screening from higher-energy states (in- 
cluding the TT bands away from the Dirac points and 
the a bands) by an additional effective static dielec- 
tric constant eint = 2.0 as done in previous graphene 
and nanotube studies [isj - lirt . The total dielectric func- 
tion e{q) is given by e{q) = 1 - v{q) [Pint(g) + P{q)] 
where v{q) = Ive^jq is the bare Coulomb interaction 
and Pyakisi) and Pici) are the static polarizabilities com- 
ing from t-h excitations involving higher-energy states 
and those involving only the low-energy tt states, respec- 
tively. Using the relation eint(g) = 1 - v{q) Pint{q) ~ Cint 
for screening with low-momentum transfer, we obtain 
e(g) « eint - v{q) P{q) ^^\^. The calculated self 
energy is then added to the LDA band energy to form 
the quasiparticle energy. Although in this scheme, the 
LDA exchange-correlation energy is not subtracted from 
the LDA band energy, it should be a reasonable approx- 



imation because the LDA exchange-correlation energy is 
nearly the same for all the tt states giving rise to a con- 
stant shift to both occupied and unoccupied states. 

We use in all the calculations a very dense grid for 
electronic state sampling corresponding to 1500 x 1500 
k-points in the irreducible wedge of the Brillouin zone 
of bilayer graphene in order to describe the extended 
wavefunction [e-h correlation length) of the excitons in 
real space, in particular at small bias voltage when the 
bandgap is small. 

The external displacement field D induces an imbal- 
ance between the charge densities on the two graphene 
layers of the BBG, which creates an internal depolar- 
ization electric field. This depolarization field induces 
additional charge changes, which in turn induce further 
adjustments in the internal electric field, and so on. We 
obtain the resulting internal electric field and the charge 
imbalance between the layers by solving Poisson's equa- 
tion ;2j. 

Figure [He) is a schematic diagram showing the 
squared e-h amplitude (wavefunction) of the lowest- 
energy optically active (bright) exciton for incident 
light with in-plane polarization, |$(re,r/i — 0)p = 

iJ^cvk^cvk (i"e|ck) (t;k|r,i = 0)|^ where the hole is fixed 
at a carbon atom belonging to the B' sublattice. The 
bound excitons [Fig. [UJa)] are comprised of interband 
transitions forming the bandgap [Fig. [ijd)]. The elec- 
tronic states in those two bands are localized at A and 
B' sublattices for the conduction and valence band, re- 
spectively, i.e., the electron and hole are localized on 
two different graphene layers. As shown schematically 
in Fig. [ije), the position of the maximum in the electron 
density is not on top of the hole, but is on a ring with a 
radius Reh from the hole that is about two orders of mag- 
nitude larger than the interlayer distance d. The radius 
of the exciton in real space is related to that in k space 
by R,h « 27r/i?c.k [see, e.g. , Fig. He) and Fig.^f)]. 

Figure [5]^ a) shows the bound exciton levels for a par- 
ticular BBG [eDd = 0.40 eV), in which we label an ex- 
citon by the radial quantum number n and the angu- 
lar momentum quantum number m of its wavefunction. 
The wavefunction of an exciton X„ „i formed from the 
free e-h pairs near the K point is approximately of the 
form <i>(re,r/i = 0) w e*'"^'''=ri™'/„_m(re) near the origin 
(re — 0) where fn,mife) has n zeros like the wavefunc- 
tions in the 2D quantum well problem having the angu- 
lar symmetry 48[. However, in our system, the angular 
symmetry is broken, i. e. , the binding energy of Xn.m 
and that of Xn,-m are different (see Table |T|. The ori- 
gin of this symmetry breaking lies in the pseudospin of 
BBG. The electronic states in in the conduction and va- 
lence bands forming the bandgap of a BBG, in the basis 
of amplitudes on the four sublattices (A, B, A' and B'), 
are 
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TABLE I: Calculated quantities of bound excitons in a BBG: 
the binding energy (-Eb), the radial quantum number n, the 
angular momentum quantum number m, and the integrated 
absorbance (lA), the absorbance integrated over energy, of 
the exciton Xn,m. made from free e-h pairs near the K point. 
The lA is for incident light with in-plane polarization. The 
quantities are the same for the exciton X'„ _„^ made from e- 
h pairs near the K' point. Here, we consider the BBG with 
Vext = eDd equal to 0.40 eV. 
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respectivelyj_where a^'s (z = 1, 2, 3, and 4) are real 



constants 



45[. As discussed above, the band edge 
states that form the bound excitons have |ai| « 1 and 
02], 1 03 1, 1 04 1 <C 1, i. e. , the electron and hole are local- 
ized at the A and B' sublattices, respectively. There- 
fore, the pseudospin of the states in a BBG imposes 
approximately an extra phase of e^^*^'' to the e-h pair 
state |ck) (wk|, resulting in an extra pseudospin angu- 
lar momentum Trips = —2. This behavior is unique in 
BBG. In pristine bilayer graphene, |ai| and |a4| are the 
same 45|, |46[ , and hence we cannot define a single extra 
phase. 

If we denote the angular momentum of an exciton 
coming from the envelope function A^^^^ by rrienv, then 
the total angular momentum quantum number (which 



is the approximate good quantum number) is given by 
m = rrienv + rn-ps- Because of the extra pseudospin an- 
gular momentum, two exciton states having Wonv of the 
same magnitude but of opposite sign are no longer de- 
generate since m would be different. Rather, two states 
having total angular momentum quantum number m and 
— m will be degenerate if the extra phase imposition by 
the pseudospin is perfect. In fact, the extra phase im- 
position of e^^'^'' is not perfect because the coefficients 
a22|, |a3|, and |a4| are non-zero, resulting in the degen- 
eracy breaking shown in Table ID [The broken angular 
symmetry shown, e. g. , in Fig. [2jm) has the same ori- 
gin.] On the contrary, due to time-reversal symmetry, 
the exciton Xn^m (formed by states near K) is degener- 
ate in binding energy with X'^ which is an exciton 
made from the free e-h pairs near the K' point with radial 
and angular momentum quantum numbers n and — m, 
respectively. Therefore, considering the spin and valley 
degeneracy and neglecting possible intervalley coupling, 
each bound exciton shown in Fig.[2][a) is four-fold degen- 
erate. 

The extra phase e"^*^'' arising from the pseudospin 
in a BBG qualitatively changes the selection rule for 
optical absorption as follows. The oscillator strength 
of an exciton S* of a BBG is given by = 
ECTk^fuk (■^'^^1 O |ck) in which O is proportional to the 
electron-photon interaction Hamiltonian. If the exciting 
photons are polarized along the x direction (i. e. , parallel 

to the graphene planes), then O cx ' '^^ ^ 



„ , where 

(7;^ 

Using Eqs. (I2|) and dS]), we ob- 
]^ (aitts e^*^'' — a2a4 e^^^). In order 



is the Pauli matrix [4] 
tain cx Eck^i , 
to have a non-vanishing oscillator strength, we should 
have oc e"^^ 



oc e 



3. , the envelope 
angular momentum quantum number TOonv should be 1 
or —1. Therefore, the total angular momentum quan- 
tum number m (which is equal to rrienv ~ 2) for the 
optically active excitons is either —1 or —3. However, 
since |ai| is by far the largest among the four |ai|'s and 
Ifliflsl ^ 10204! [4^, effectively, only the excitons Xn,-i 
or X'^ I are optically active (Table I]). This unusual op- 
tical selection rule in a BBG, hence, originates from the 
unique pseudospin physics. 

In the discussion below on the optical absorbance, for 
concreteness, we shall assume that the polarization of 
the incident light is linear and is parallel to the graphene 
planes. Accordingly, the lowest-energy exciton Xq.o 
[Figs. [2][b)-[2][d)] is dark and the second lowest-energy 
exciton Xo._i [Figs. [2je)-[2jg)] is bright. As seen from 
the calculated oscillator strength in Table HI the lowest- 
energy bright excitons by far dominate the absorbance 
spectrum. The first, second and third bright excitons 
have zero, one, and two nodes in the exciton wavefunc- 
tion along the radial direction, respectively, in both mo- 
mentum and real space (Fig. [J). Also, there are many 
dark exciton levels between the bright exciton ones as 
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shown in Fig. [UJa). A change in the polarization direc- 
tion of the incident hght away from the graphene plane 
would alter the optical strength of the levels from those 
given in Fig. HJa). 

In a 2D hydrogen atom, the binding energy is pro- 
portional to (n + \m\ + l/2)~^ resulting in a 2N + 1-fold 
degeneracy with iV = n -I- |m| As shown in Table [H 
however, this degeneracy in the binding energy of the 
excitons in a BBG is broken, and, further, the order of 
the binding energies largely deviates from the case for a 
2D hydrogenic model. Also, we have checked that the 
detailed order of exciton levels changes with the external 
displacement field. 

Figure [3] shows the calculated absorbance spectrum of 
BBG (for in-plane linearly polarized incident light) near 
the bandgap energy and the wavefunction of the lowest- 
energy bright exciton that forms the main peak for sev- 
eral bias voltages. Remarkably, when e-h interactions 
are accounted for, the absorbance profile is dominated 
by a single four-fold degenerate excitonic level with huge 
oscillator strength. Accordingly, the dominant feature 
of the absorbance profile near the bandgap energy be- 
comes symmetric when excitonic effects are considered 
- as in carbon nanotubes [13, E^; whereas, if these ef- 
fects are neglected, highly asymmetric absorbance spec- 
tra are obtained reflecting the "effective" ID van Hove 
singularity in the joint density of quasiparticle states dis- 
cussed above. The huge excitonic effects observed here in 
fact originate from this ID singularity [s^l which becomes 
more and more dominant as the bandgap increases. On 
the contrary, excitonic effects on the low-energy (< 1 eV) 
optical response of pristine bilayer graphene are negligi- 
ble since its joint density of states is characteristic of a 
2D system '30|. The enhancement of excitonic effects 
with the bandgap is reflected in the increase in the ex- 
citon binding energy [Fig. IDJa)] and the decrease in the 
exciton radius [Fig. [3] and Fig. IDJb)] . 

In a previous study (4ll | , we have shown that when the 
photo-excitation energy is close to the energy of the zone- 
center optical phonons in BBG (~0.2 eV), Fano line- 
shapes in the absorbance profile develop due to the cou- 
pling of e-h pair excitations with the phonons. We expect 
that similar exciton-phonon coupling behavior, whose ef- 
fects on the optical response is large when the optical en- 
ergy gap is around 0.2 eV, would arise if electron-phonon 
interactions are taken into account. 

The above results are applicable to suspended 
BBG 50]. However, for BBG on substrates, excitonic 
effects are altered due to enhanced screening from the 
substrate. As an example, we consider the effect of back- 
ground screening due to the substrate on the optical re- 
sponse of BBG relevant for the experimental setup in 
Refs. |2^ and ^ij. For substrates above and below the 
BBG having dielectric constants ei and £2, respectively, 
their effect can effectively be replaced by a single mate- 
rial having a dielectric constant of eec = (ei + e2)/2 [5l[ . 



Using the static dielectric constant of Si02 (=3.9) and 
that of amorphous AI2O3 (=7.5), we may roughly set 
the external background dielectric screening as eec = 
(3.9 -I- 7.5)/2 = 5.7. Figure [5] shows similarly calculated 
quantities as in Fig. lUJa), but now for BBG with added 
substrate screening as discussed above. The exciton bind- 
ing energy = A*^^ — A^^^ is smaller than the case 
without substrate screening. The calculated optical gap 
^BSE jg reasonable agreement with the experiment 
A'^'^P [23. 

In this work, the inter-layer trigonal warping effects on 
the electronic structure (owing to an atom on one layer 
interacting with further neighbors on the other layer), 
i. e. , the trigonal anisotropy in the energy band disper- 
sion near a Dirac point ^4], have been neglected. If these 
effects were taken into account, the cylindral symmetry 
about an individual Dirac point is weakly broken, lead- 
ing to only minor changes in exciton energies and to some 
of the dark excitons gaining very small optical oscillator 
strength. However, the change in the overall absorbance 
spectra (which are dominated by excitons) at different 
gaps is negligible [H^l- 

In conclusion, we have shown that excitons in bi- 
ased bilayer graphene dramatically change the optical 
response because of the ID nature of the joint density 
of quasiparticle states in this system. These excitonic 
effects are remarkably tunable by the external electric 
field. Also, we have shown that the pseudospin character 
of the electronic states dramatically alters the excitonic 
structure (energy level degeneracy, optical selection rule, 
etc.) of this system. These results illustrate the richness 
in the photophysics of biased bilayer graphene and their 
promise for potential applications in nanoelectronic and 
nanophotonic devices at room temperature 27 1. 
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FIG. 1: (a) Schematic diagram showing the structure of pristine bilayer graphene whose unit cell is composed of four different 
sublattices (A, B, A', and B'). (b) Schematic bandstructure of pristine bilayer graphene (origin is the Dirac point). Solid blue 
and dashed red lines represent valence bands and conduction bands, respectively, (c) and (d): Same schematic diagrams as 
in (a) and (b) for bilayer graphene under a displacement field D generated through a double-gate. In (d), A is the energy 
bandgap and vertical arrows represent interband transitions responsible for the formation of excitons. (e) Schematic diagram 
showing the probability density that a photo-excited electron is found at when the hole (blue empty circle) is fixed at the 
origin, |'&(re, — 0)^ (see text). For visualization purposes, we show the quantities in a vertical plane that includes the hole. 
The fake thickness of the plotted profile (red) is proportional to the probability density. The interlayer distance d is extremely 
exaggerated in (e). The size of the exciton R^h is much larger than d [Fig. |3{b)]. 



FIG. 2; (a) Calculated free e-h pair excitation dispersion (-Bck — Suk versus k) and exciton levels of a BBG with an external 
electrostatic potential between the two graphene layers Vcxt = eDd equal to 0.40 eV (Fig. [T]). Thick red lines and thin blue 
lines show optically active (bright) and inactive (dark) exciton levels, respectively, for incident light with in-plane polarization. 
The exciton Xn,m (-'^n.-m) formed by e-h pairs near the K (K') point is denoted by its radial quantum number n, angular 
momentum quantum number m (see text), and binding energy i5b. Each exciton level is four-fold degenerate due to the spin 
and valley degeneracy (see text). There are many other higher-energy bound excitons not shown here whose energy is below the 
bandgap. (b) The squared amplitude of the lowest-energy exciton [exciton Xo,o in (a)] in momentum space |^f„k| . (c) Squared 
wavefunction in real space of the corresponding exciton in (b). The plotted quantity is the probability density |<l"(re, — 0)^ 
of finding an electron at Ve given that the hole is fixed at one of the carbon atoms (at the center of the figure) in sublattice 
B' (Fig. [T]). (d) Real part of the exciton wavefunction Re $(re,rh = 0) for the corresponding exciton in (b). (e)-(g), (h)-(j), 
and (k)-(m): Similar quantities as in (b)-(d) for the first, the second and the third bright excitons [excitons Xo,-i, and 
X2,-i in (a), respectively]. 



FIG. 3: (a) Calculated absorbance spectra of BBG (with an arbitrary energy broadening of 5 meV and in-plane polarization) 
where Kxt = eDd (see Fig. [l]) is 0.14 eV. Results with (blue or solid line) and without (red or dashed line) e-h interaction 
effects are shown, (b) Wavefunction of the lowest-energy bright exciton (Jfo,-i or Xq i) that forms the dominant peak in the 
absorbance spectrum. The plotted quantity is the probability density \^{re,rh = 0)|^ of finding an electron at re given that 
the hole is fixed at one of the carbon atoms (at the center of the figure) in sublattice B' (see Fig. [1]). (c) and (d), (e) and (f), 
and (g) and (h): Same quantities as in (a) and (b) for Vcxt = 0.27 eV, 0.40 eV, and 0.66 eV, respectively. 



FIG. 4: (a) The quasiparticle bandgap A*^^, the optical bandgap A'^^'^, and the binding energy (= A'^'' - A^^'^) of BBG 
versus Kxt = eDd. (b) The size Reh, defined in Fig. [T][e), of the lowest-energy bright exciton (Xo.-i or Xq i) versus Kxt. The 
line is a guide to the eye. 



FIG. 5: The quasiparticle bandgap A*^^, the optical bandgap A'^^'^, and the binding energy Et, (= A^^^ - A'^^'^) of BBG 
under background screening (ebg ~ 5.7) versus Kxt ~ eDd. Measured data A^^'' are taken from Ref. [S^j. 



